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Abstract 

Operator  formalism  is  employed  for  introduction  of  effective  constitutive  tensors  of  bian- 
isotropic multilayered  periodic  strucures  in  a wide  wave  band.  It  is  based  on  the  approxi- 
mate calculation  of  the  characteristic  matrix  of  the  unit  cell  of  the  system  with  the  help  of 
Campbell-Hausdorff  series  [1].  In  this  paper  we  show  that  effective  constitutive  tensors  can 
be  correctly  introdused  for  a greate  variety  of  mediums. 


1.  Introduction 

One  of  the  most  fruitful  approaches  in  the  theory  of  propagation  of  waves  in  multilayered  media  is 
the  operator  formalism  based  on  the  use  of  intrinsic  representation  of  vectors  and  operators  (see 
[2],  [3],  for  example).  It  enables  one  to  find  compact  coordinate-free  formulae  for  very  complex 
systems  and  to  avoid  the  cumbersome  calculations  required  by  usual  components  techniques  and 
caused  mainly  by  the  necessity  of  adopting  a different  coordinate  system  for  each  different  layer. 

In  this  paper  we  use  this  formalism  to  find  the  effective  constitutive  tensors  of  the  medium 
formed  by  a periodic  set  of  plane  bianisotropic  layers  with  different  thicknesses  ln  (n  = 1, 2, . . . N, 
where  N is  the  number  of  the  layers  which  constitute  the  unit  cell). 

In  what  follows  we  assume  the  constitutive  relations 

Dn  = £n^n  4"  anHjii  Bn  = /3nEn  + Hn&n,  (1) 

for  layer’s  materials,  where  £n,fin  and  an,  are  the  dielectric  permittivity,  the  magnetic  per- 
meability tensors  and  the  pseudotensors  of  gyrotropy,  respectively.  These  equations  can  be  also 
written  in  matrix  form,  as 


Dn 


(2) 


If  e„,  /in,  an,  f3n  (n  — 1, 2, . . . N)  are  complex  nonsymmetric  tensors,  then  equations  (1)  or  (2) 
describe  an  absorbing  anisotropic  and  gyrotropic  medium,  subject  to  the  influence  of  external 
electric  and  magnetic  fields  and  elastic  deformations. 

The  main  purpose  of  this  paper  is  to  introduce  the  medium’s  effective  material  tensors, 
feasible  for  the  use  in  the  wide  wave  band. 
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2.  Effective  Constitutive  Tensors  in  Wide  Wave  Band 


To  treat  this  problem  we  use  here  an  approach  to  introduce  the  effective  material  parameters, 
proposed  in  [4].  It  is  based  on  the  approximate  calculation  of  the  characteristic  matrix  P of  the 
unit  cell  of  the  system  with  the  help  of  Campbell-Hausdorff  series  [1]. 

The  characteristic  matrix  P = exp  ( ikolM ) of  a layer  with  thickness  l relates  the  six-vectors 
(E,H)t  at  the  layer  boundaries  (fco  = cj/c).  The  matrix  M can  be  written  in  the  form  [5] 


M = rqQx  ( n-Bx)Tq , 


(3) 


Mo"  ”tf  ) ’ B*  = ( b»  T)’ 


T„  = £-Nq,  S 


, A rq  = K~lQ{n-Bx),  Q = 


0 ) 

q®q  / 


(4) 


(5) 


where  operator  TZq  1 is  the  inverse  to  the  operator 


K„=  ( £gl  ^ 

9 ^ 0g  1 fiql  ) ’ 


(6) 


eq  = qeq,  Hq  = q/«q,  aq  = qaq,  j3q  = q/?q, 

qx  is  the  antisymmetric  dyadic  dual  to  the  unit  normal  of  the  boundaries  q,  b is  the  tangential 
component  of  refraction  vector  m [2], [3],  q ® q is  the  dyad  and  1 is  the  unit  dyadic. 

In  the  case  of  a medium  formed  by  two  alternate  layers  with  different  thicknesses  ln  and 
different  sets  of  tensor  constants  en,pn,o;n,/3„  (n  = 1,2),  the  characteristic  matrix  of  the  unit 
cell  has  the  form 

P = exp  ( ikoLM ) = exp  exp  (zfc0ZiMi) , (7) 

where  L = + 1%  is  the  system  period  and  M is  some  matrix.  This  matrix  can  be  expressed  in 

terms  of  the  layers  parameters  with  the  help  of  Campbell-Hausdorff  series  as  follows 


M 


■d’h 


fiMy  + f2M2  + i/,^  [M2,  Mi]  - 


(8) 


where  fn  — ln/L  is  the  relative  thickness  of  the  n-th  layer,  fn  = 1,  [M2,My]  = M2Mi  — 

M\M2.  If  Trln/\  <C  1,  then  the  series  (8)  quickly  converges  and  one  can  drop  the  remainder  of  it 
after  the  fc-th  term.  Here  for  the  sake  of  simplicity  we  have  limited  our  consideration  by  three 
terms  of  series.  The  first  and  the  second  terms  correspond  to  the  long  wavelength  limit. 

The  tensor  parameters  A,B,C,D  can  be  directly  used  for  the  analysis  of  the  wave  propa- 
gation in  the  considered  system,  for  example  for  finding  the  reflection  and  transmission  tensors 
[3]. 

For  the  sake  of  simplicity  we  assume  q to  be  the  left  and  the  right  eigenvector  of  each  of 
the  tensors  en,  fin,  an,  /3n  (enq  = qen  = enqq,  • • •) . Besides,  we  shall  use  in  (7)-(8)  the  matrices 
M/,  Mu,  M2/,  describing  the  transformation  of  the  tangential  components  of  the  field  vectors, 
instead  of  M, Mi,M2,  describing  the  transformation  of  the  full  six- vectors.  In  the  majority  of 
applications  this  would  be  ample.  Under  above-mentioned  conditions: 


a <g>  a 

0 


0 

a <g>  a 


r, 


Mni  = Qx 


(9) 
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r= 


I = 1 - q ® q, 


where  TZ^q  is  the  transposed  operator,  a = bqx , 


[M2I,Mu}  = Qx[R21  - A(bl)]T, 


(10) 


(11) 


722i  = ^Qx'R-i  — 72iQx722, 


(12) 


A (a)  = Al 


a ® a 0 
0 a ® a 


+ 


a ® a 0 
0 a ® a 


(13) 


where 


Al 


72f0 


7?t 

KuQx^1 


S2 


> An  = 


721,  72T0 

^<3x721/  -~^-Qxn2I 
0 2 <?1 


(14) 


and  72n/  = T72nT,  n = 1, 2. 

For  the  correct  introduction  of  permittivity,  permeability  and  gyrotropy  tensors  of  the 
medium  equivalent  in  the  considered  approximation  to  the  system  at  hand,  the  matrix  Mj 
is  bound  to  depend  on  the  effective  tensors,  angle  and  plane  of  incidence  in  the  same  way  as  in 
the  homogeneous  medium.  From  (8)— (13)  it  follows  that  these  conditions  are  met  if  A(a)  = 0. 
In  this  case  we  can  find  the  effective  material  tensors  in  the  same  way  as  in  [6] 


" ) = r (/i72i  + f2K2  + i/i^7221)  r + 7 lgQ, 


(15) 


with 

Kq^ihK^  + h^y1.  (16) 

For  example,  A(a)  = 0 if  A l = A#  = 0.  That  takes  place  provided  [6] 


722/72-29  = 72i/72i9,  TZ2q7l2 / = 72ig72i/, 


(17) 


Note,  that  both  relations  (17)  are  equivalent  if  ai  — /%  = 0.  Even  so  it  is  rigid  enough  require- 
ment, signifying  that  tensors  £u  and  £2/,  pu  and  p2 / bound  to  have  the  same  eigenvectors. 

Of  course,  above  mentioned  conditions  need  not  be  met  exactly  having  regard  to  approximate 
nature  of  method.  We  can  find  valuations  ||Al||  and  ||722i||  instead.  It  should  be  more  than 
sufficient  if  ||Al||,  1 1 A /?|  | far  less  than  ||722i || - In  what  follows  we  shall  use  Euclidean  valuation 


||X||  = [(XXt)t]5,  (18) 

where  is  Hermitian  transposed  tensor  and  (XX^)t  is  a trace  of  tensor  XX^ . 

By  way  of  example  let  us  consider  nongyrotropic,  nonmagnetic  layers  with  real  symmetric 
tensors  £*.  In  this  instance 


IIAill2  = 


2 

+ 2sin  4? 


(19) 


||722i||2  — 2 [(£2+  — £i+)2  + (^2-  — £i-)2  + 2sin <j>2  (ei+  — £i-)  (£2+  — £2-)]  , (20) 

where  4 is  the  angle  between  the  preferred  axes  of  crystals,  a+,  £j_  (i  = 1, 2)  are  the  eigenvalues 
of  tensors  eu  = l£il-  Formulae  for  magnetic  crystals  are  similarly  to  (19), (20).  Since  usually 


134 


|| a* ||,  HAH  <C  ||£ill)  IlMill)  corresponding  formulas  are  valid  also  for  analysis  of  the  gyrotropic 
layers. 

Inasmuch  a2  ~ e,9,  from  (19),  (20),  (13)  it  follows  that  for  identical  layers  (£i+  = £2+, £i-  = 
£2-)  values  ||A(a)||  and  ||^2i||  are  approximately  equal.  Situation  changes  in  the  case  of 
different  parameters  of  the  layers.  In  this  case  it  is  possible  to  choose  parameters  of  the  crystals 
providing  the  fulfillment  of  condition 


II A (a)  II  « 11^21  II-  (21) 

Relation  (21)  can  be  provided  much  easier  for  structures  comprising  more  then  three  layers  in 
the  unit  cell.  Then 

M = f\M\  + /2M2  + /3-M3  + (22) 

ITT 

j {hh  [MSi  M2]  + fsh  [Ms,  Mi]  + hh  [M2,  Mx}}  + . . . , 

and  instead  of  ||  A (a)  ||,  ||7?-2i||  (12)— (14)  we  shall  have  more  complicated  but  more  versatile 
expressions,  which  can  be  changed  by  varying  thicknesses  of  the  layers.  It  should  be  noted  that 
sometimes  the  ’’wide  wave  band  approximation”  reduces  to  the  long  wavelength  approximation. 
In  particular,  if  TZ\  = 7Z3,  l\=h  = h,  then  [M3,  Mx]  = 0,  [M3,  M2]  = - [M2,  M{\  and  structure 
do  not  possess  properties  of  form  gyrotropy  [6]. 

Formulas  (15),  (9)-(12)  are  valid  for  arbitrary  material  parameters  of  the  layers  for  a small 
enough  angles  of  incidence  and,  of  course,  in  the  case  of  normal  incidence  ( a = 0). 

3.  Conclusion 

It  is  shown  that  effective  constitutive  tensors,  not  depending  on  the  angle  and  plane  of  incidence 
(true  constitutive  tensors),  can  be  introduced  in  the  wide  wave  band  for  the  great  variety  of 
anisotropic  periodic  mediums.  Formulas,  which  enable  one  to  determine  whether  it  is  possible  to 
introduce  true  constitutive  tensors  for  the  structure  at  hand  or  not,  are  obtained.  In  the  frame 
of  considered  approximation  it  is  hardly  probable  to  introduce  correctly  constitutive  tensors  for 
an  arbitrary  parameters  of  the  layers.  From  the  other  hand,  at  normal  incidence  there  are  no 
limitations  on  the  parameters  of  the  layers  and  obtained  effective  tensors  may  be  very  useful  for 
synthesis  of  new  materials. 
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